双参数混合正则化方法及在电导率反演成像中的应用          下载免费PDF全文

王文娟  Chris Farmer  John Ockendon  冯宝宾  曹俊兴  谭永基 《地球物理学报》2011,54(8):2154-2159

通过引入带有二阶正则算子的正则化项,建立了一种双参数混合正则化方法.为确定最佳正则化参数,这里主要应用L-曲线法、偏差原理和广义交叉校验准则的优化组合来确定.首先对理论模型进行了数值模拟,通过与截断奇异值分解法、共轭梯度法及标准Tikhonov正则化法的结果比较,表明该方法不仅精度高,而且对于数据的随机扰动具有稳定性....  相似文献   

An algorithm for fast 3D inversion of surface electrical resistivity tomography data: application on imaging buried antiquities     

N.G. Papadopoulos  P. Tsourlos  C. Papazachos  G.N. Tsokas  A. Sarris  J.H. Kim 《Geophysical Prospecting》2011,59(3):557-575

In this work a new algorithm for the fast and efficient 3D inversion of conventional 2D surface electrical resistivity tomography lines is presented. The proposed approach lies on the assumption that for every surface measurement there is a large number of 3D parameters with very small absolute Jacobian matrix values, which can be excluded in advance from the Jacobian matrix calculation, as they do not contribute significant information in the inversion procedure. A sensitivity analysis for both homogeneous and inhomogeneous earth models showed that each measurement has a specific region of influence, which can be limited to parameters in a critical rectangular prism volume. Application of the proposed algorithm accelerated almost three times the Jacobian (sensitivity) matrix calculation for the data sets tested in this work. Moreover, application of the least squares regression iterative inversion technique, resulted in a new 3D resistivity inversion algorithm more than 2.7 times faster and with computer memory requirements less than half compared to the original algorithm. The efficiency and accuracy of the algorithm was verified using synthetic models representing typical archaeological structures, as well as field data collected from two archaeological sites in Greece, employing different electrode configurations. The applicability of the presented approach is demonstrated for archaeological investigations and the basic idea of the proposed algorithm can be easily extended for the inversion of other geophysical data.  相似文献   

2010年1月12日海地 Mw7.0地震InSAR同震形变观测及同震滑动分布反演          下载免费PDF全文

薛莲  孙建宝  沈正康 《地震地质》2011,33(1):157-174

2010年1月12日GMT时间21时53分,在海地境内(72.57°W,18.44°N)发生了Mw7.0地震.文中利用干涉合成孔径雷达(InSAR)方法获得了覆盖整个震区的高精度形变观测资料,用以研究该地震的发震机理.采用ALOs PALSAR数据,分析了轨道、大气等误差源对干涉信号的影响,最终获得了雷达视线向(LOS...  相似文献   

Determination of gravity anomaly at sea level from inversion of satellite gravity gradiometric data   总被引：1，自引：0，他引：1  

Mehdi Eshagh  Lars E. Sjberg 《Journal of Geodynamics》2011,51(5):366-377

Gravity gradients can be used to determine the local gravity field of the Earth. This paper investigates downward continuation of all elements of the disturbing gravitational tensor at satellite level using the second-order partial derivatives of the extended Stokes formula in the local-north oriented frame to determine the gravity anomaly at sea level. It considers the inversion of each gradient separately as well as their joint inversion. Numerical studies show that the gradients T_zz, T_xx, T_yy and T_xz have similar capability of being continued downward to sea level in the presence of white noise, while the gradient T_yz is considerably worse than the others. The bias-corrected joint inversion process shows the possibility of recovering the gravity anomaly with 1 mGal accuracy. Variance component estimation is also tested to update the observation weights in the joint inversion.  相似文献   

Application of constrained sparse spike inversion in reservoir predication:a case study of Oriente Basin in South America     

WANG Qing  LU Zhanguo . Key Laboratory of Geo-detection of Ministry of Education  China University of Geosciences  Beijing  China 《东北亚地学研究》2011,(2)

In the conditions of low Signal-to-Noise Ratio(SNR) of seismic data and a small quality of log information,the consequences of seismic interpretation through the impedance inversion of seismic data could be more precise. Constrained sparse spike inversion(CSSI) has advantage in oil and gas reservoir predication because it does not rely on the original model. By analyzing the specific algorithm of CSSI,the accuracy of inversion is controlled. Oriente Basin in South America has the low amplitude in geological...  相似文献   

基于细化分层法探讨面波频散曲线反演参数的简化     

邵广周  李庆春 《地球科学与环境学报》2011,(3):317-320

常规频散曲线反演过程中需要不停地改变分层数、层厚度和层速度等参数,实现过程相对繁琐,而采用细化分层法对反演参数进行简化则避免了上述缺点。具体思路为:根据目的层探测深度(如20m)将地下介质分为若干个(20个)厚度为1m的薄层和1个均匀半空间层(共21层),这样在反演中分层数和层厚度均为已知参量,反演过程只需修改速度参数即可,避免了改变分层数和层厚度等参数,显著简化了反演计算过程。正演计算和反演结果均表明:细化分层与实际分层计算出的频散曲线是等效的,细化分层反演结果的总体效果与真实模型非常接近,这说明细化分层方法用于频散曲线反演是切实可行且有效的;将地下介质划分为1m厚的薄层,反演后每层均可得到1个横波速度,能满足反演分辨率的要求;由于实际地下介质的速度是随深度渐变的,细化分层后比按频散曲线拐点分层(每分层的厚度可能是几米或几十米,同一分层内介质的横波速度相等)更接近实际情况。  相似文献   

塔河油田3区石炭系超深薄层碎屑岩储层预测研究   总被引：2，自引：0，他引：2  

盛海波  刘蕊  付国民  蔡玥  张庭娇 《地球科学与环境学报》2011,33(4):390-396

针对塔河油田3区石炭系卡拉沙依组地震数据很难识别薄砂体的问题,对比分析了采用地震资料提频处理、分频处理、储层高分辨率敏感参数反演等不同手段得到的储层预测结果,筛选出高分辨率储层敏感参数反演作为卡拉沙依组砂泥岩段储层主要预测手段.通过探讨砂泥岩对声波时差的响应,以及对比分析自然电位、自然伽马以及补偿中子孔隙度与声波时差的关系,确定自然电位与声波时差的相关性最好,由此选取自然电位作为最佳电性敏感参数参与高分辨率敏感参数反演预测.从井点、剖面、平面等方面检验和评价高分辨率储层预测成果,结果表明其反演预测效果较好,符合塔河油田3区石炭系的储层预测要求,实现了对石炭系主力砂体空间展布的解释与描述,建立了塔河油田3区石炭系卡拉沙依组储层模型.  相似文献   

民勤大气边界层特征与沙尘天气的气候学关系研究   总被引：4，自引：0，他引：4  

李岩瑛  张强  薛新玲  王荣基 《中国沙漠》2011,30(3):757-764

 为了更好地理解西北干旱区大气混合层(ML)厚度的变化特征及其对当地沙尘气候形成的影响, 利用民勤2006—2008年3—6月逐日08时和20时探空资料、降水和日最高气温,计算和分析了最大混合层厚度、逆温层特征和垂直风场及其对沙尘气候形成的影响。结果表明,民勤沙尘天气的大气边界层有显著的昼夜变化,白天厚、逆温强而多;沙尘天气的最大混合层厚度在2 600 m左右,介于无降水与有降水天气之间;扬沙主要由锋面中的冷空气引起,而沙尘暴主要由低层风场的剧烈扰动和500 m以上高层冷锋入侵引起。沙尘暴发生前近地面风场有明显的扰动,沙尘暴发生时在500 hPa以下有显著的冷空气活动,白天较强。能见度小于100 m的强沙尘暴夜间风速大,冷空气较强。  相似文献   

层序地层学在隐蔽圈闭预测中的应用--以松辽盆地南部为例   总被引：8，自引：0，他引：8  

柏涛  徐志伟 《吉林大学学报(地球科学版)》2004,34(1):73-78

根据松辽盆地南部地区典型的反演地震剖面及典型井的体系域、层序的发育特征,在松辽盆地南部地区共识别出9个层序、12个体系域,并应用层序地层学理论对松辽盆地南部地区青山口组及姚家组隐蔽圈闭进行了预测.研究结果表明:砂岩透镜体岩性圈闭主要分布于青二段及姚二三段高水位体系域;砂岩上倾尖灭岩性圈闭主要分布于姚一段水进及高水位体系域;地层超覆圈闭主要分布于各个水进体系域;断层-岩性圈闭和构造-岩性圈闭可以分布于各个体系域,分别取决于断层和构造发育层位.  相似文献   

Palaeoproterozoic mid‐basin inversion in the northern Mt Isa terrane,Queensland     

P. G. Betts 《Australian Journal of Earth Sciences》2013,60(5):735-748

The Mellish Park Syncline is located in the northern part of the Mt Isa terrane. It has an axial trace that transects the remnants of the unconformity‐bounded Palaeoproterozoic Leichhardt and Isa Superbasins. The syncline is separated into a lower and upper component based upon variation in fold geometry across the basin‐bounding unconformity. The lower syncline, in the Leichhardt Superbasin, is tight and has an inclined west‐dipping axial plane. The upper syncline, in the Isa Superbasin, is open and upright. The geometry of the lower syncline is a consequence of a period of shortening and basin inversion which post‐dated the Leichhardt Rift Event (ca 1780–1740 Ma) and pre‐dated the Mt Isa Rift Event (ca 1710–1655 Ma), forming an open and upright north‐oriented syncline. Subsequent southeast tilting and half‐graben development during the Mt Isa Rift Event resulted in the lower syncline being tilted into its inclined geometry. Sequences of the Isa Superbasin were then deposited onto the eroded syncline. The geometry of the upper syncline reflects regional east‐west shortening during the Isan Orogeny (ca 1590–1500 Ma). The position of the upper syncline was largely controlled by the pre‐existing lower syncline. At this time the lower syncline was reactivated and tightened by flexural slip folding.  相似文献